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The temperature dependence of paramagnetic susceptibility of the superfluid He-B in aerogel 
is found. Calculations have been performed for an arbitrary phase shift of s-wave scattering in 
the framework of BCS weak coupling theory and the simplest model of aerogel as an aggregate 
of homogeneously distributed ordinary impurities. Both limiting cases of the Born and unitary 
scattering can be easily obtained from the general result. The existence of gapless superfluidity 
starting at the critical impurity concentration depending on the value of the scattering phase has 
been demonstrated. While larger than in the bulk liquid the calculated susceptibility of the B- 
phase in aerogel proves to be conspicuously smaller than that determined experimentally in the 
high pressure region. 

PACS numbers: 67.57.Pq, 67.57.Bc 



Superfluid '^He is an ideal object for a study of physical 
properties of a system with non-trivial pairing. In par- 
ticular, it is interesting to investigate how a superfluid 
state could be influenced by the presence of impurities. 
Direct contamination of '^He with any atomic impurities 
is impossible. Last years instead of this experimental- 
ists use Hquid '^He to flll up aerogel which is a matrix of 
randomly arranged silica fllaments of nanometer diame- 
ter. Commonly used aerogels occupy about 2% of space 
volume, the rest being taken with liquid helium. 

Experiments ||i|, |||, H, ||||, |||| estabHshed decrease 
of the superfluid transition temperature and the density 
of the superfluid component of ^He in aerogel. These 
effects are related to the scattering of quasiparticles 
from surface, which suppresses Cooper pairing in the p- 
state. Corresponding theoretical treatments have been 
proposed in papers H , R] , §1 , ||9| , - 

NMR measurements fj pointed on the existence of the 
A-phase of ■^He in aerogel in relatively high magnetic 
flelds as well as of the B-phase [^ in lower fields. Re- 
cently, evidence of the phase transition between the A- 
and B-phases has been demonstrated iQ, [|3|, ||l^ . 

The A-phase as an equal spin pairing state possesses 
practically the same spin susceptibility as the normal 
phase. This property remains intact in the A-phase in 
aerogel. On the other hand, in the B-phase where all 
the three spin states of the Cooper pairs with the spin 
projections Sz = 0, ±1 are equally populated, spin sus- 
ceptibility is partly suppressed compared to its normal 
state value. This suppression in aerogel was experimen- 
tally found |[l2J to constitute about 45% of the value in 
the bulk 3He-B. 

The purpose of the present article is to calculate spin 
susceptibility of the '^He-B in aerogel. We shall work in 
the framework of BCS weak coupling theory, aerogel be- 
ing approximated as homogeneously distributed ordinary 
impurities. At high pressures this model gives values of 
A and ps roughly by a factor of two larger |q| , [i| than 



those observed experimentally. At the same time, as it 
was discussed in paper ||l§|, at low pressures, when the 
coherence length is larger than the average pore size, the 
homogeneous scattering model is more likely to be justi- 
fied. 

A similar conclusion follows from the present deriva- 
tion. The calculated susceptibility is less suppressed com- 
pared to the bulk dependence. This suppression is less 
pronounced for larger values of the scattering phase but 
still even in the unitary limit noticeably larger than the 
measured suppression of the susceptibility at high pres- 
sures ll, |l2[. 

The calculations presented in the paper have been done 
for an arbitrary value Sq of the phase shift of s-wave 
scattering. Both extremes: the Born approximation and 
the unitary limit can be easily obtained from the general 
result. 

Corresponding theory of changes of paramagnetic sus- 
ceptibility in ordinary superconductors caused by impu- 
rities with ordinary, paramagnetic and spin-orbital type 
of scattering has been developed in the Born approxima- 
tion in paper ||l^. The present theory is not a simple 
generaHzation of the cited paper. The point is that usu- 
ally to come out of the Born approximation in the theory 
of dirty alloys it is sufficient to substitute the Born scat- 
tering amplitude by its exact counterpart [|^ . This sub- 
stitution does not lead to the correct answer for a dirty 
metal in a magnetic field. The method developed in the 
present paper does not rely on the substitution. 

The paper is organized as follows. After the introduc- 
tion of the basics of the Abrikosov-Gor'kov theory in the 
first section, we write down its solution for an arbitrary 
impurity scattering phase without magnetic field in the 
second section, where we also estabHsh the region of the 
existence of gapless superfiuidity. Then by finding cor- 
rections in leading order in magnetic field we arrive at an 
expression for the susceptibility. In the last section we 
briefiy discuss conclusions. 



I. EQUATIONS 

The Abrikosov-Gor'kov equations of an electrically 
neutral impure superfluid in an external magnetic field 
B are (see, e.g. |p^) 

(lUJn - H{k) - SK)) G(k, u;„) = i, (1) 

where w„ — ■KT[2n + 1) is the Matsubara frequency. 
Green function G'(k, a;„) is a 2x2 matrix in particle-hole 
space 

G(k,a;„) - (^ ^+^(k,c.„) G„,(k,c.„) ) ' ^^^ 

made up of the normal Gaf3{'k,0Jn) and anomalous 
i^Q^(k, [jjn) Green functions which are in their turns 2x2 
matrices in spin space. Here 

G„/3(k,t^„) = -G;^^(-k,-a;„) (3) 

and the superscript T implies transposition. 
The Hamiltonian 



( i?0a/3(k) A^ 



(4) 



'ka(3 "Ual3\ 

consists of the one-particle part 

^Oa/3(k) = ^k<5a/3 ^ M'^a/jB (5) 

including kinetic 

a = Clkl = A:V2m* - ep (6) 

and Zeeman energy, and of the pairing interaction 
Vk/3,A/^(k, k') via the order parameter A^,^^, set by the 
self-consistency equation 

^kcp - -TY,J2^f3c.M^,k')F,xiu;^,k'). (7) 

n k' 

B is the magnetic flux, (Tap is a vector of the three Pauli 
matrices in spin space. 

The impurity scattering self-energy part comprises nor- 
mal Si and anomalous S2 parts 



E = 



Si S2 

Sn Si / ' 



(8) 



and in the clean limit of small impurity concentrations 
n-imp, when one can neglect the interference of the scat- 
tering on different impurities, is just rUmp times the con- 
tribution of a single impurity 



^ ^iinp-' 1 



(9) 



where T is the scattering amplitude. It is related to 
the scattering potential C/kk' by the Lippmann-Schwinger 
equation 



Tkk' — C^kk' + /_^ C^kqGqTqk' 



(10) 



In the simplest case of a short-range impurity poten- 
tial U{r) = u6{r) its Fourier component is momentum- 
independent C/k-k' — u, so that 



U, 



kk' 



Uk-k'Sai3 







-c/k'-k5.J="'^"^^^' (11 



where r is a vector of the three Pauli matrices in particle- 
hole space. Parameterized in terms of the scattering 
phase So the impurity potential is 



u = — taiiSo/TTNo 



(12) 



Eq. (^ is then easily solved, and the scattering am- 
plitude turns out to be independent of momenta 



fkk'('^n) = T{uJn) = ull- C/^G(k,C^„) j 



(13) 



Equations (0), (||), (||), ( [T^ ) form a closed system of 
equations on the Green function (g). 

The pairing interaction Vj3a,\^ii}^,'k') is usually reck- 
oned non-zero only in a thin ~ ei <C ef layer near the 
Fermi surface. In "^He spin-orbital interaction is weak and 
the Cooper pairs are formed in the spin-triplet, orbital 
p-wave state. So the pairing interaction can be factorized 



Vaa,AA<(k,k') = --Fifcfc'g/3agJ^, 



(14) 



where Vi is the constant of the p-wave pairing attraction 
and 



gQ/3 = iicr(Ty)ai3 

is a vector of the three basis symmetric matrices. 
The order parameter 



A, 



kaf3 



dfcga/3, 



(15) 



(16) 



where d ? = A^iki is the order-parameter vector. 
Substituting (|l|) and ^) into (§) yields 

dfe = ^V,TY,Y.Ckk')i^As^F{cu,„-k% (17) 



n k' 



As one can see 



d,d* = -tr.[AtA,], 



(18) 



where trcr denotes a trace over spin indices. For unitary 
phases, among which are the A- and B-phases actually 
being realized in pure ^He, {AfA~)ap c>c 6ap and so it 

is convenient to introduce the notation A- = d-d-. We 

k k k 

will consider only the unitary phases. 
In the A-phase 



(19) 



d^. -^/^Ad[(A' + zA")fc], 



where d, A', and A" are three unit vectors, A' and A" 
are mutually orthogonal. 
In the B-phase 



dj^ = ARke"^, 



(20) 



where R is the matrix of rotation on an arbitrary angle. 

In both cases, the scalar A is introduced so as to fulfill 

the normalization (A-)- — (d-d!)- — A^, where {■ ■ ■): 

stands for the averaging over directions k of momentum. 



II. ZERO MAGNETIC FIELD 

In the absence of magnetic field the solution of Eqs. 
(0), (|), (|), (111) is (see elsewhere) 



g(0) 



1 



cJi, 



^k + A| 



(21) 



-{iuJn + S,k)SaP 



kap 



At 

fea/3 



'{iuJn -^k)^Q/3 



The corresponding sum over momenta in the expres- 
sion ( [l^ ) for the scattering amplitude is 

^G(")(k,u;„)=.g(^„)W„/3, 

k 

where we denoted 

k 

= —inNouj. 



Here Nq = m*kF/2TT^h'^ is the density of states on the 
Fermi surface. 

So the scattering self-energy M) equals 




£(0) 



''imp •- 



% + ugl 



0af3- 



(23) 



1 — u'^g'^ 

In the above formulae both the Matsubara frequency 
LUn ~ 7TT{2n+ 1) and the kinetic energy are renormalized 

iun = it^„ + -tr^[S(°)(tj„)], (24) 

(25) 



1 



a = ek-^tr,[?3S(")(c^„) 



Here trr is the trace in particle-hole space. 

It should be noted that the term proportional to ra 
in E = £(•') 4- S(^) + ... must be omitted. This term 
just produces a shift in the chemical potential (renor- 
maHzation of ^k ) from introducing the impurities and 
disappears in the assumption of particle-hole symmetry 
after integration over ^ in the subsequent calculations. 



Substituting (H) and (H) into (H) yields the self- 
consistency equation on the renormalized Matsubara fre- 
quency UJn 



LUr. 



^l + ^l 



^^l + cos^6oAl/^ 



where T is the scattering rate 



nNo 



■ 2 c 

sm Oq. 



(26) 



(27) 



In the Born limit {Sq -^ 0) it tends to the conventional 
half of the inverse free flight time 



'imp r2 



TrA^n 



6o = nin^pTiNou- 



1 

27' 



(28) 



In the B-phase A? = A^ does not depend on k, one 

can thus omit angular brackets of the fc-averaging in (E6h . 
Using the overt form (|2l|), the self-consistency Eq. p^ ) 
taking the trace over spin indices and averaging over k 
reduces to 



1 



NoVi 



-^E 



d^ 



Lot, 



A' 



--^E 



1 



(29) 



This equation determines the temperature and impu- 
rity concentration behavior of the order parameter A. To 
obtain an expression of value one should exclude the un- 
observable pairing constant Vi off the left-hand side of 
Eq. (pg]). Using the standard procedure [l^, one obtains 
the self-consistency equation on A(r, F) 



1 ^ 



ttT 






V^ 



—- . (30) 



The asymptotics of A(T, F) at T 
Appendix. 



0, Tc are found in the 



Expanding (Eq) in degrees of A gives w„ 



Fsgncjn + O(A^). When inserted into (30) it yields in 
zeroth order in A^ the Abrikosov-Gor'kov ||l^ implicit 
expression for the dependence of the critical temperature 
Tc on the impurity concentration 



1 "^cO , 

log — — w , 

^ Tr V 2TTTr 



^-a 



(31) 



where i/" is the digamma function. Note however that F 
corresponds to different impurity concentrations for var- 
ious scattering phases (see (|2^)). T'c(F) vanishes at the 
critical concentration 



Tc = \A 



00 J 



where 



iOO 



TrTcoe 



(32) 



(33) 
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FIG. 1: Order parameter A versus temperature T for dif- 
ferent scattering phases So and impurity concentrations ob- 
tained numerically. Solid lines correspond to the Born limit 
So = 0, dashed lines to So — j and dotted lines to the unitary 
scattering So = ^. In order of the decrease of the superfluid 
transition temperature Tc the triples of the curves correspond 
to r/Fc =0, 0.2, 0.4, 0.6, 0.75, 0.9, 0.98. The first curve in 
the pure case F = is common for all phases. 



is the order parameter at zero temperature for a pure 
superfluid, 7 ~ 0.5772 is Euler's constant. 

Fig. |l| shows temperature dependences of the order 
parameter for 7 impurity concentrations and for three 
scattering phases: Sq = 0,j, and ^ obtained numerically 
from (|3C|). Non-shown plots for intermediate < Jq < f 
all sit in the interior of the belt between curves for So ~ 
(Born limit) and (5o = f (unitary limit). 

The order parameter Aq at zero temperature as a func- 
tion of the scattering rate F/Fc is presented on Fig. g 
for several values of the scattering phase 60 ■ The curves 
for intermediate values of < (5o < f lie inside the strap 
between the curves for the Born {Sq -^ 0) and unitary 
{Sq -^ ^) limits. 

From the Abrikosov-Gor'kov theory of paramagnetic 
impurities in alloys it is known that in the Born limit 
{So -^ 0) starting at the scattering rate of 



2r, 



^7r/4 



o.9ir^ 



(34) 



the energy spectrum gap vanishes though the amplitude 
A of the order parameter remains flnite until Tc. Ordi- 
nary impurities in a p-wave superfluid will give rise to the 
existence of an analogous gapless region of impurity con- 
centrations. Indeed, density of states of quasiparticles at 
an energy e off the Fermi level is expressed in terms of 



the Green function as 

7V(e) = --35]G(k,c.„)|,^„_ 

TT ^ — ' 
k 

Using ( p2| ) we get 

N{e) = iVo3 



»e+iO ■ 



VA2 - ^ ' 



(35) 



(36) 



where e — iujn\i^ -^e+io- From ( pq) we obtain an implicit 
equation on e for a given e: 



e = e + r 



?VA2 - ?2 

A^ COS^ Sq —'P 



(37) 



For e = it gives three solutions in the complex plane: 
£1=0 and 



2A2 cos2 So-T^ - F Vr2+4A2^ii?^ 



2,3 



(38) 



It can be verified that for all the three solutions A^ — e^ 



is always (real) positive, and so a/A^ 
order for N{e 



- always real. In 
0) to be finite one than sees from (136 



that e should be negative. From ( |38[) this gives 
F > Acos^(5o. 



(39) 



This inequality means that for a given phase shift Sq 
and a scattering rate F the superfiuid is gapless at 
temperatures Tgi < T < Tc, where Tgi is such that 
F > A(Tgi) cos^ (5o. The superfiuid is gapless in the 
whole temperature range < T < Tc oi the existence 
of superfiuidity if inequality ( p9| ) is fulfilled already for 
Ao = A(T = 0). 



Making use of the implicit dependence (A7) of Aq on 
Sq, we get for the lower boundary of the onset of gapless 
superfiuidity in the whole temperature range 



F = 2F^ cos^ ^1 



^^^ Sq 



oe 



(40) 



which is a generalization of the Abrikosov-Gor'kov value 
(Q) on an arbitrary scattering phase ^o- The correspond- 
ing value of the upper boundary of the order parameter 
at zero temperature is 



Aq = AooP/ 



(41) 



When the scattering phase is increased (impurity po- 
tential becomes "stronger") the boundary value of F de- 
creases, meaning that starting from less impurity concen- 
tration the superfiuidity becomes gapless. At last, in the 
unitary limit Sq ^ ^ the boundary value is F = 0, and 
the liquid is gapless for whatever small impurity concen- 
tration. 

The points of the onset of gapless superfiuidity are 
marked on the plot of the dependence of the order pa- 
rameter at zero temperature Aq/Aqo on the suppression 
Tc/Tco of the superfiuid transition temperature (see Fig. 
g) . The regions of the existence of gapless superfiuidity 
are shown on the plots of Ao(T'c) with bold lines. 
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FIG. 2: Order parameter Ao at zero temperature as a func- 
tion of the suppression Tc/TcO of the critical temperature for 
several values of the scattering phase Sq. For each So the value 
of Tc/Tco is pointed, the suppression below which makes the 
superfluid gapless in the whole range < T < Tc. In the 
unitary limit {do ^ f) the superfluid is gapless for whatever 
Tc/Tco- Bold parts on each curve mark the regions of the 
existence of gapless superfluidity. 



III. MAGNETIC SUSCEPTIBILITY 

Magnetization M can be calculated from the normal 
Green function in a magnetic field B as 



M = flT^^Gal3(Tpa. 



(42) 



n k 



To find the static susceptibility Yq/j (Mq, = Xa/3-B/3) it 
is sufficient to solve Eqs. (0), (Q), (g), jl^ ) in leading 
order in the field. Let the superscript (l}designate the 
quantities proportional to B. From (||), (|3|) we get 



g(i) = ^s(o)c/(^gW)c/-is 



(0) 



(43) 



One should plug here the expression for the Green 
function in first order in the field derived from (fu) 



GW=G(°)(ifW+EW)G(°). 



(44) 



We are now going to show that both the anomalous 
self-energy E2 remains zero and the order parameter A 
unchanged in first order in magnetic field: 



.(1) 



0, 



A(i)=0. 



(45) 



:a lly, this IS 
solution of (p^), (14^ ) produces (Hnear) homogeneous 

equations on S2 ^^'^ ^^^^ which have regularly only 
trivial solutions. At the same time, the equation on the 
normal part of self-energy 'Ei involves a right-hand side 
proportional to the magnetic field coming from 



H'^^Hk) = 



(1) + 
'kaP 



A 



(1) 

kaf) 



(46) 



To be more specific, when written out in components, 
Eq. (H) reads 



G(^) - G(") [1] + i^(°) [2] , 



G<^2], 



(47) 
(48) 



where we made use of the abbreviated notations 

[1] = (E(^)-Ai<TB)G(") + (4'^-A(i))F(")+, (49) 
[2] = (E(i)+-A(i)+)G(")H-(El''+/.<TTB)F(o)t50) 

Let us first look on the equation on A*-^'. Expanding 
Eq. (0) yields 



A^^'+ = ;Fig+r5:^(fcfc')ti.[ff,F(i)+(c.„,k')]. 



n k' 



(51) 



Substituting here ( p8| ) for F^^^+ one obtains the equa- 
tion on A'-^^^. It turns out that after taking trace over 
spin indices as well as integrating over ^ on the supposi- 
tion of electron-hole symmetry there remain only terms 
proportional to A*^^)"*". 

Indeed, consider the terms in (Q) in pairs. The two 

terms involving the anomalous part Sj ('^n) of the self- 
energy, which, we remind, does not depend on momenta, 

are proportional to either G G^-^\ which does not de- 
pend on k, or to (f (0)+)^ ex k^ . They vanish after aver- 
aging with k' over the direction of momenta in (O) . 

The two terms proportional to B when inserted into 
(|l|) give 

ex tr,[5^ (f(")+<tBG(") - G^"Vtbj^(o)+^] 

ex d^ tr^[5^ {gtii^n + £,')cr\ - crj(i(l;„ - C')5i^)]-Ba- 

The terms proportional to ^' vanish on integrating over ^' 
in the supposition of particle-hole symmetry. The terms 
proportional to a;„ vanish on summation over the Mat- 
subara frequencies w„ = (2n -t- l)7rT from n = —00 to 00 
because of the oddness of a;„(wn): 2„(— Wn) = — a;„(w„). 
Similar argumentation leads to the conclusion of van- 
ishing of the terms involving Y^l because, as we shall see 
later, in the assumption of particle-hole symmetry E*^ 
is even in w„. 



So we are left with a homogeneous linear equation on 
A*-^-* which has normally only trivial solutions A^^' = 0. 

One obtains the equations to determine the normal 
T,\ and anomalous S2 parts of the self-energy by plug- 



ging (21) and (M) into 



^(1) - 



.(1) 



-(1) 



E«+ S^ 



^inipf- 



(1 - uV) 



2„2\2 



(52) 



"V)Ek^('^ 






Since A^ , and consequently also i^*-"' (a;„, k) by virtue 



.(1) 



of (21), and A^ are all proportional to the unit vector 



k, any of the terms in ([17[), ( |4q ) comprising one or three 
of either of the quantities vanishes on summing over k 
in (p3|). For the anomalous part there remains only two 
terms giving a homogeneous equation on the four com- 
ponents of S2 : 



-(1)+ 



y2„2 A^ 
^ k 



/'^(o)+j.(i)^(o)+ _^^(o)j.a)+(^(o) 

(53) 



Unless in exceptional cases, this equation has only trivial 
solution E2 — 0- However, one can verify separately for 
each component of the expansion of S2 (a;„) in terms of 
the three Pauli matrices and the unit matrix: 



.(1 



.(1) 



(1)/ 



^2q/3('^") ~ ^2 {^n)Saf3 + ^2 {^n)f^al3 



(54) 



that this homogeneous equation on Sj {ujn) does indeed 
have only trivial solution. 

Taking into account the remark after Eq. (|2^) we omit 
the terms proportional to T3 in (p2) and for the normal 
part Til get the equation 



Si 



(1) 



^impW 



2 2 
U 9 



(i-"V)^^ 



[G^"\t[^' ~ H(tB)G 



(0) 



^(o)^v^i) 



+ F(^'{T\''+ficT'B)F'-">+] 



(55) 



Makin g us e of the explicit form of the zero-order Green 
function (|2l|), we first calculate the sum over k, introduc- 
ing the auxiliary quantity 



n(c^„) 



TTiVnA^ 



3((:D2 + A2)3/2 



(56) 



k,w -; - ^nK) and ^^F(")F(")- 



so that Ek(G*°^)' = I 
n(cj„). In addition, we introduce 



n(a;„) 



1 



n 



u^g^ 



imp •- 

rA2 



(1 - uV)^ 
Zof. cos 2(5o 



n(c^„) 

+ A^ cos^ 5o 



3 ^SJTA^(w2+A2cOs2,5o) 



2 • 



(57) 



The second terms in ( pq ) is proportional to the prod- 
uct of Pf} ex Agaa^B^k and F('')+ cx Ag+g^k. After 



averaging over fc and integration over ^ we find that 

^F(")/(c.„)F(«)+ = in(a;„)g/(c.„)g+. (58) 



So 

sw = ^n(c.„) 



^r 



3~/ 



S^^' - mtB + ig(sl'^ + MT^B)gt59,) 
Sj'^ + a*<tTb + ig+(s(^) - /irTB)^0,) 



Expanding e![ and S^ in the basis of the three Pauli 



matrices and the unit matrix similarly to (|5J), one sees 
that for the coefficient of expansion over the unit matrix 
one obtains an uncoupled homogeneous equation, with 
trivial solution. So one may put simply 



.(1) 



•^1 cr, 



-(1) 



Using g — icray, g+ = —iaycr and the equalities 



ay era,. 



we find that 



So that solution of the system (p9[), (pfl) is 



SmCtSm = dt^^g^i^ = o-^- 



(61) 
(62) 
(63) 



Sl^) = -;.<tB- "K) 



1 - n(c^„) 



s1^^=m-tb "K) 



Since from (|47|) it follows that 



1 - n{ujn) 

(64) 



AirTTB)i^(o)+, 



we find from (^2f) that 

M 



^ 1-n 



(Wn) 



(65) 



(66) 



Because we have first integrated over k and left the 
summation over n this expression does not include ||lj 
normal state Fermi gas PauH susceptibility x^ ^ 2^'^Nq. 
Adding and subtracting it from both sides of ( |4^ ) , we get 
eventually 



A2 



Xl 3 ^ 



(^2 + A2)3/2 



n 1 



rA2 



52 cos 2(5o + A2 C0s2 (Sq 



3 V^n + A2 (52 + A2 C0S2 ,5o)' 

(67) 



This expression gives susceptibility as a function of tem- 
perature for a given scattering rate T and a scattering 
phase (So- The infinite sum over the Matsubara frequen- 
cies should be taken using self-consistent expressions for 
the renormalized Matsubara frequency cj„ (J26|) and the 
order parameter A(T) (|3y). 

In the clean limit F ^- the Matsubara frequencies rest 
unrenormalized and susceptibility reduces to its conven- 
tional bulk value 
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(69) 



is the B-phase Yosida function. 

At r — > Tc susceptibility x"/Xn behaves linearly. In- 
deed, then A — > and 2 — > tj + F. Hence 
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(70) 



where Xc = F/27rTc and ^c"^ = V''"^! + ^^c) is the 
polygamma function. Using the asymptotic expression 
( |A2| ) for A at r — > Tc, we find eventually 
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In the Born limit {5q -^ 0) we arrive at the expression 
for the spin susceptibility 
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(72) 



(73) 



Susceptibilities for several impurity concentrations and 
three scattering phases ((5o = — Born limit, 5q — ^, 
and (5o = f — unitary limit) are plotted as functions of 
the reduced temperature T/T^ on Fig. ^ with the Fermi 
liquid effects taken into account (see below). 
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FIG. 3: Spin susceptibility of superfluid He in aerogel versus 

^ temperature for several values of impurity concentration F 

and scattering phases (5o with the Fermi liquid corrections 

taken into account. Solid lines correspond to the Born limit 

(71) ^0 = Oi dashed lines to So — ^ and dotted lines to the unitary 



scattering 5o 



In order of the decrease of the superfluid 



transition temperature Tc the triples of the curves correspond 
to F/Fc =0, 0.2, 0.4, 0.6, 0.75, 0.9, 0.98. The first curve in 
the pure case F = is common for all phases. 



A. Fermi-liquid corrections 

Fermi-liquid interaction between quasiparticles leads 
to renormalization of the susceptibility. An external mag- 
netic field is screened by a polarization of the liquid, 
which is quantitatively described by even terms of the 
expansion in Legendre polynomials of the antisymmetric 
(exchange) part F^ of the Fermi-liquid interaction. 

It was argued in Ref. [ pO[ on the basis of experimental 
data that only the zeroth term F^ is non-zero. Then the 
magnetization Ma = XapBp induced by magnetic field 
B in a liquid with interaction equals the magnetization 
M, 
free liquid by an effective field 



X^flSfl^ that would be induced in an interaction- 



B<=ff = B - F^M/xl 



(74) 



Here Xn = 2/i^iVo is the Pauli paramagnetic susceptibil- 



Aq/3- 



ity of a normal Fermi-gas. This yields ||2J, 

In the B-phase x° o ex 6ai3 and we get 
X _ l + F§ x° 



Xr: 



1 I pa X" yO ' 
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(75) 



(76) 



where Xn = Xn/(l+^(f ) is the Fermi-liquid susceptibility. 



In ^He Fn^ 



slightly varying with pressure. 



For the asymptotics at T ^ Tc we have 



X_ 

Xn T~^T( 

where from d7i 



dix/Xn) 



d{T/T,) 



T-Tr 



T^T^ 



Tr 



d{x/Xn) 



diT/T,) 



1 



rf(x7x^) 



T^Ta 



l + F^ diT/T, 



(77) 



(78) 



T-^T^ 



is the slope of the plots of x versus T at T ^ T,. 

Taking — | as a value for Fq , we see that the slope of 
the curve x (T) at T — Tc with the Fermi-Hquid interac- 
tion taken into account, is 



1/(1 + ^o) 



(79) 



times greater than that of x" (T) — without the interac- 
tion. 

The slopes of both x° (T) and x (T) for F^ = -^ ver- 
sus the suppression Tc/Tcq of the critical temperature are 
plotted on Fig. |j in two scales of the ordinate axis — the 
left one for the slope of x° (T) , and the right one for the 
slope of X (T), which differs from the former in ( [79|) times. 
On the same picture we also pointed the slopes experi- 
mentally observed in ||^ at 18.7 bars and in [Q at 32 
bars — the only experimental data present up to day the 
authors are aware of. The error bars span the possible 
slopes of the Hnear fits to the data at T <Tc within the 
plotted error bars of the respective experimental works. 



IV. CONCLUSIONS 

In the framework of the homogeneous scattering model 
we found spin susceptibility of ^He in aerogel considering 
the latter as ordinary randomly distributed impurities 
with an arbitrary scattering phase So in s-wave channel. 
The answer is given in the form of an infinite sum over the 
Matsubara frequencies w„ = (2n+l)TrT with the summa- 
tion term depending on the amplitude A of the order pa- 
rameter and the self-consistently renormaHzed Matsub- 
ara frequency cj„. Numerically plotted curves x(F) show 
that susceptibility is less suppressed in comparison to the 
bulk dependence, this suppression being less pronounced 
for larger Sq. 

Comparison with high-pressure experimental data 
(18.7 bars |§ and 32 bars ^) shows that the slope of 
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FIG. 4: The slope at T ~» Tc of the reduces temperature 
T/Tc plots of the "bare" susceptibility x /Xn (l^ft scale) or the 
susceptibility x/Xn with the Fermi liquid corrections (right 
scale, coefficient of proportionality (JTOh) as functions of the 
suppression Tc/Tco of the critical temperature. The crosses 
mark the experimental values of the slopes in |B| (D) and in 

@(0). 



x(r) at r -^ Tc is approximately twice less than the 
value predicted by the theory. The discrepancy is obvi- 
ously unavoidable because of the major inadequacy of the 
homogeneous scattering model at high pressures, where 
the superfluid coherence length becomes shorter than the 
correlation length of the internal structure of aerogel ||l5[ . 
Note that quantitatively the difference between experi- 
mental and theoretical values seems to be in Hne with 
the fact that experimentally observed at high pressures 
(19.7 bars) superfluid density |^^ was approximately by 
the factor of 2 smaller than predicted by an Abrikosov- 
Gor'kov theory [g. 

On the other hand, we pointed out that for any fi- 
nite scattering phase So > the impurity concentration 
threshold of the gapless superfluidity onset diminishes 
dramatically. For example, for the value So = j which 
was maintained to be the most satisfactory phase for 
the real structure of aerogel, gapless superfluidity sets 
on in the whole temperature range starting at the crit- 
ical temperature suppression of T, ~ 0.51Tco compared 
to Tc « 0.22Tco in the Born limit Sq -^ 0. Whereas in the 
unitary Hmit Sq ^ ^ the superfluid is always gapless. For 
impurity concentrations exceeding the gapless threshold 
all the thermodynamical quantities vanish algebraically 
at zero temperature. Experimental observation of such 
behavior could provide additional immediate insight into 
the quasiparticle scattering on impurities in aerogel. 
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where the divergent integration up to infinity is cut off 
at ei. 

The general temperature dependence of the order pa- 
rameter is obtained by combination of (P9|), (33) and ( [Aq ) 



1 



1 



log-^ = 7rTy (■ ^ ^ 
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APPENDIX A: ASYMPTOTICAL BEHAVIOR OF 
THE ORDER PARAMETER A AT T ^ Tc AND 

In this Appendix we elaborate on the details of the 
calculation of the asymptotics of the order parameter al- 
ready utilized in the main text of the paper. 



1. Asymptotics at the critical temperature 



At r ^ the value of the order parameter tends to 
the value Aq the implicit dependence of which on F is set 
by replacing summation in the above expression with an 
integral. The i nteg ral over lu from the first term in the 
parentheses in ( [Aq ) transforms into an integral over uj by 
means of the substitution (|^ and eventually we get an 
equation that determines implicitly the impurity concen- 
tration behavior of the value of the order parameter at 
zero temperature: 



Expanding ( |2q ) in the degrees of A^ gives (we con- 
sider only positive w„ because of the oddness of ojni^n)'- 

ii„ « u;„ + r + rl-^^^A^ + o(a4). (ai) 



When inserted into ( po| ) it yields in first order in A^ 
the asymptotics at T ^ Tc 
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Here Xc = T/2t:Tc and tpr' = ^^"^5 + ^c), where 



V(") {z) = d-Mz) = Y. T^^I^ (A3) 



fc=0 



[z + k) 



is the polygamma function. At F ^ also Xc -^ and 



,,(2) 



^If-i ^ V^^Hi) = -14C(3), and thus (^) transforms 
into the conventional BCS asymptotics 
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if F > Ao cos^ 5. Here 



^o(r,A) = 



' F2 - 2 A2 cos2 5o + F Vf2+4A2^h?5^ 

(A9) 



is the solution of ( 26 ) for LOn = and in ( |A8D it was taken 
with the arguments F, Ao. 



For small F Eq. (A7) gives 



2. Order parameter at zero temperature 

Consider the sum ttT Y^ni^l + A^)-!/^. When T ^ 
the summation over the Matsubara frequencies w„ may 
be substituted with an integration over the continuous 
variable uj 



^-^ . /rxTL A2 
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while for F close to the critical value 
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